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Using high-throughput Tethered Particle Motion single molecule experiments, the double-stranded 
DNA persistence length, Lp, is measured in solutions with Na"*" and Mg^"*" ions of various ionic 
strengths, I. Several theoretical equations for Lp{I) are fitted to the experimental data, but no 
decisive theory is found which fits all the Lp values for the two ion valencies. Properly extracted 
from the particle trajectory using simulations, Lp varies from 30 nm to 55 nm, and is compared 
to previous experimental results. For the Na"*" only case, Lp is an increasing concave function of 
7“^, well fitted by Manning’s electrostatic stretching approach, but not by classical Odjik-Skolnick- 
Fixman theories with or without counter-ion condensation. With added Mg^+ ions, Lp shows a 
marked decrease at low 7, interpreted as an ion-ion correlation effect, with an almost linear law in 
7“^, fitted by a proposed variational approach. 

I. INTRODUCTION 

Ions play a major role in the cell, for example by mod¬ 
ifying the protein activity, inducing a voltage between 
the intra-cellular and extra-cellular matrix, or control¬ 
ling the DNA packaging in viral capsids or in the nu¬ 
cleus. Complex mechanisms involving DNA such as its 
wrapping around histones or its denaturation will only 
be thoroughly understood when the role of mobile ions 
on the DNA conformation gets elucidated, DNA being 
one of the most charged biopolymers found in nature. 

The first quantitative experimental studies of DNA 
conformational properties as a function of salt concen¬ 
tration have been done in 1978 by Harrington [T], us¬ 
ing flow birefringence (FB) experiments, to measure the 
DNA radius of gyration in dilute DNA solutions. The 
DNA radius of gyration is intimately related to the DNA 
persistence length, namely the correlation length of the 
tangent-tangent correlation function 

(t(s) • t(0)) = exp {-s/Lp) (1) 

where t(s) is the unit vector tangent to the chain at 
the point of curvilinear index s. The persistence length, 

Lp, thus characterizes the chain stiffness at small length 
scales. 

Experimentally, the persistence length cannot be di¬ 
rectly measured and the required procedure to extract 
it has remained a major issue since these first quanti¬ 
tative measurements [T]. Other optical techniques have 
then been used such as the transient electric birefringence 
(TEB) [2] or magnetic birefringence (MB) [^, linear 
dichroism (LD) [5], dynamic light scattering (DLS) [3], 
and force-stretching by optical tweezers (EOT) [7] were 


used to estimate the variation of Lp as a function of the 
ionic strength 7. In a recent paper, Savelyev [1] reviewed 
the available experimental data and showed that they 
could be divided into two groups, based on the distinct 
behaviours of Lp found at high ionic strength. Indeed 
whereas the first group of experimental data [UliHT] indi¬ 
cated a slow decrease of Lp with increasing 7, the second 
one [21HHE] found a significant one. 

Hence no global picture emerges yet from the litera¬ 
ture. Many reasons can be put forward, such as the diffi¬ 
culty to estimate accurately the ionic strength in buffers, 
which is not simply equal to the added salt concentra¬ 
tion, or the method of extraction of the persistence length 
from the experimental observables. Indeed, the extrac¬ 
tion of the variations of Lp with 7 by using FB, MB and 
LD techniques is very sensitive to the optical arrange¬ 
ment between the electric field and the molecular axis, 
which is related to the tedious evaluation of the mag¬ 
netic/optical anisotropy of a single base-pair. Moreover, 
in LD experiments, the mechanism of interaction between 
nucleic acids and the electric field depends on the polar¬ 
izability of the ionic cloud surrounding the DNA, which 
therefore requires an additional modeling. In FOT ex¬ 
periments, the high DNA stretching induced by the force 
modifies the DNA structure and the organization of its 
ionic cloud. Finally in DLS experiments, the DNA hy¬ 
drodynamic radius is estimated from the measure of the 
diffusion coefficient. Inferring the persistence length is 
not easy especially due to excluded volume effects. 

Among the different theories developed to explain the 
variations of Lp with the ionic strength, 7, the most fa¬ 
mous is the Odjik-Skolnick-Fixman (OSF) [131 [14] theory 
where Lp is the sum of a bare, non-electrostatic persis¬ 
tence length and an electrostatic contribution scaling as 
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. This theory accounts, at least qualitatively, for the 
fact that a rise in I provokes an increasing screening of 
the repulsive phosphate ions of the backbone which leads 
to a more flexible DNA chain. Taking into account the 
so-called Manning counter-ion condensation around the 
DNA at low J, the prefactor in front of , propor¬ 
tional to the square of the effective DNA charge, is low¬ 
ered. All these experimental works have therefore been 
compared to these types of approaches, with no deci¬ 
sive conclusion [4]. More recently, Manning proposed 
a new theory [T^, and explicitly considers the electro¬ 
static stretching force of the polyelectrolyte, which qual¬ 
itatively fits some experiments [7] or numerical results at 

high I mile]. 

In this paper we reconsider the old issue, still under 
debate of the dependence of DNA conformation, at 
room temperature, with the ionic strength of the sur¬ 
rounding solution, using the recent single-molecule tech¬ 
nique of high-throughput Tethered Particle Motion (HT- 
TPM) min]. We measure the persistence length of two 
DNA of 1201 and 2060 base-pairs (bp) for a large range of 
well controlled ionic strengths from I = 10“^ to 3 mol/L 
with Na’*' counter-ions and with or without Mg^+ added 
counter-ions. In the first Section, we present the sim¬ 
ple and well controlled HT-TPM experiments and their 
analysis. The next Section is devoted to the experimental 
results and the extraction of the persistence length from 
the HT-TPM amplitude of motion using numerical simu¬ 
lations. Our results are then compared to previous ones, 
and to the various existing theories as well as our de¬ 
tailed variational approach and an interpolation formula 
that fits all our experimental Lp. Finally our concluding 
remarks are given in the last Section. 


II. EXPERIMENTAL SECTION 

A. High-Throughput Tethered Particle Motion 
(HT-TPM) experimental procedure 

DNA molecules were produced by PCR amplifica¬ 
tion (oligonucleotides from Sigma-Aldrich) : Biot-F1201 
5’-CTGGTGAGTACTCAACCAAG-3’ and Dig-R1201 
5’-CTAGAATGGATGGGAAGG-3’ on pTOGl plasmid 
and Biot-F2060 5’-GTGGAATGATAGGGGGAGAC-3’ 
and Dig-R2060 5’-TGAGTTGGGCGTTTGGAGAC-3’ 
on pBR322. 

HT-TPM permits the simultaneous tracking of hun¬ 
dreds of single DNA molecules free to fluctuate in so¬ 
lution that are tethered to a coverslip at one end and 
labelled by a 300 nm particle at the other end (see 
Fig. 0. HT-TPM on chip assembly is performed as pre¬ 
viously described in detail in Ref. (TH]. The anchoring 
of the DNA-particle complexes to the neutravidin (Invit- 
rogen) printed sites was performed in PBS buffer (Eu- 
romedex) supplemented with Img/mL of pluronic 127 
(Sigma-Aldrich) and 0.1 mg/mL BSA (Sigma-Aldrich), 
noted T-BSA-Plu. 



FIG. 1: Principle of the Tethered Particle Motion experi¬ 
ment. The DNA is tethered at one end on the coverslip and 
labelled by a particle (radius Rp) at the other end. The pro¬ 
jected particle position R|| is tracked as a function of time. 
The DNA is model by a chain of beads in the simulations. 


A large range of buffers were then used to test the effect 
of both ion valency, by using monovalent Na^ or divalent 
Mg^"*", and salt concentration, on the DNA conformations 
(Tables 1 and 2 in the Supporting Information). The 
first buffer, corresponding to zero salt added, and named 
Zero-salt-buffer in the following, is a phosphate buffer 
(KH 2 PO 4 1 mM, Na 2 HP 04 3 mM, pH 7.4, pluronic F127 
1 mg/mL). From this one, we added successively various 
concentrations of NaCl or MgCl 2 to obtain a large range 
of salt conditions (X-salt-buffer). Before starting the ex¬ 
periment, the flow cell was extensively rinsed with the 
Zero-salt-buffer (~ 100 chamber volumes), let incubated 
during 1 h at room temperature, then rinsed again with 
~ 100 chamber volumes of Zero-salt-buffer. The experi¬ 
ment started with a Zero-salt-buffer measurement, next 
the concentration in monovalent ions was progressively 
increased by addition of ^ 100 chamber volumes of X- 
salt-buffer. Then the flow cell was extensively rinsed with 
the Zero-salt-buffer (~ 100 chamber volumes), incubated 
during 4 min and rinsed again with by ~ 100 chamber 
volumes with the Zero-salt-buffer. A new Zero-salt-buffer 
measurements was performed. At last, the divalent ion 
concentration was progressively increased and new acqui¬ 
sitions were carried out. For all conditions, the acquisi¬ 
tions were performed at a controlled temperature equal to 
25°C, and a movie of 1 min were recorded and analyzed. 
We ensure the reliability of the experimental procedure 
by checking the agreement between the two values of the 
Zero-salt measurement obtained before the addition of 
monovalent ions and before the addition of divalent ions. 
Experiments were repeated on different days to ensure 
the reproducibility of our results. 

The tethered particles of diameter 300 nm were visual¬ 
ized using a dark-field microscope (Axiovert 200, Zeiss) 
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equipped with a x32 objective and an additional xl.6 
magnification lens and a temperature control system 
(Physitemp TS-4MPER). Images were acquired during 
1 min at a frame recording rate of 25 Hz on a CMOS 
camera Dalsa Falcon 1.4M100. The field of observation 
covers an area of 215 /im x 160 /im. 


B. HT-TPM procedure of analysis 

The software developed by Magellium (France) tracks 
in real time the positions of all the particles, corrects for 
experimental drift, calculates the asymmetry factor to 
select tethered particles valid for the analysis, and finally 
the experimental root mean square end-to-end distance 
projected on the surface, i?oxp|| = the amplitude 

of motion of the bead, along the time trace (see Fig. [^. 
We invite the reader to refer to Ref. m for the detailed 
calculations of i?exp||- 

In order to measure the small differences expected on 
Rcxplp we set up a two-step procedure for the analysis of 
i?cxp|| where, first, a criteria of validity, then some correc¬ 
tions are applied. All this procedure was performed with 
homebuilt Mathematica scripts. Details can be found in 
Ref. [T5]. During this procedure, around 12% were elim¬ 
inated and the final number of kept trajectories for each 
DNA condition typically ranges between 100 and 1000 
(see Tables 1 and 2 of the Supporting Information). 
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III. EXPERIMENTAL RESULTS 

A. Simple estimate of DNA persistence length 
from the amplitude of motion 


FIG. 2: Experimental HT-TPM amplitude of motion vs. 
the ionic strength I for a DNA of length (a) L = 2060 bp and 
(b) L = 1201 bp with Na”*" (red) and with added Mg^”*" (blue) 
cations. Vertical lines delimit the close-to-physiological salt 
conditions, as listed in Table [I] 


In Figures and b are shown the variations of the 
experimental amplitude of motion, i.e. the root mean 
square end-to-end distance projected on the surface, 
Rcxpip measured by HT-TPM, as a function of the ionic 
strength, I, of the various buffers listed in Tables 1 and 2 
(Supporting Information) for two DNA of lengths 2060 
and 1201 bp. The ionic strength is defined as 

( 2 ) 

i 

where Zi and Ci are respectively the valency and the con¬ 
centration of ion i in the buffer. The red (resp. blue) 
symbols correspond to the monovalent Na"*" ions (resp. 
with divalent Mg^+ ions added) buffers. The black ones 
correspond to Zero-salt-buffer conditions. 

The plots clearly show that, when I is raised, i?oxp|| 
decreases of about 10%, for monovalent as well as diva¬ 
lent ions. This is thus a subtle effect. This behavior is 
qualitatively due to the fact that negatively charges car¬ 
ried all along the DNA, repel themselves more strongly, 
when the ionic strength is low, this repulsion decreas¬ 
ing when the ionic strength increases due to electrostatic 


screening. This self-repulsion increases the rigidity of the 
macromolecule thus increasing i?cxp||- Furthermore, one 
notices that the variation of i?exp|| is amplified when the 
DNA length, L, is increased. This is consistent with the 
fact that the process under examination occurs homoge¬ 
neously along the DNA. 

Comparing the experimental results for monovalent 
ions and divalent ones on the same I range, one observes 
that the decrease of Roxp||(^) is more important with 
Mg^+ ions, once again an effect amplified for long DNAs. 
Moreover, one notes a sudden drop of i?oxp|| at very low 
concentration of Mg^+ (from I = 10.1 mmol/L, without 
Mg^+, to 10.5 mmol/L). However, it is difficult to know 
with high precision the ionic strength for such low val¬ 
ues, since we cannot exclude that some additional ions 
get released from the surface, for instance. Therefore, 
the errors bars are potentially large for these points. 

To access to properties of the DNA itself, it is necessary 
to correct the effect of the particle and the glass substrate 
to access to properties of DNA only, such as its end-to- 
end distance, Rdna = \/ (R-dna)- simplest way is to 

substract the particle radius by assuming that the DNA 
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extremity and the particle, of radius Rp = 150 nm, move 
independently, which leads to [19] 

(RDNA) = ^(R||)-i?p (3) 

By doing so, we neglect the effect of the glass substrate. 
It has been taken into account analytically only for long 
and flexible polymers |23|, which is not the case for the 
DNA of this study. Indeed their length {L = 2060 and 
1201 bp, i.e. 700 nm and 408 nm, using 1 bp= 0.34 nm) 
is no more than a few persistence lengths {Lp « 50 nm), 
which allows us to qualify them as semi-flexible. As al¬ 
ready shown in [11] the approach proposed in [53] does 
not work well for such DNA lengths [66] . 

In Figures la and b of the Supporting Information are 
shown the DNA end-to-end distance, i?DNA, obtained 
by Eq. ([^, as a function of the inverse of the ionic 
strength, I~^. Of course the relative effect of the salt 
is slightly higher, once the particle radius is deduced. 
The next step consists in extracting the DNA persistence 
length Lp from IJdna- The simplest way would be to use 
the Worm-Like Chain (WLC) formula, valid for a phan¬ 
tom chain in solution, resulting from Eq. Q [53] : 

(R?)NA)=2L^(;^-l + e-^/^^) (4) 

However, this way of extracting Lp leads to quite high 
values of the persistence length compared to the com¬ 
monly accepted values around 50 nm. For instance at 
low I, Lp saturates around 76 nm for the 2060 bp long 
DNA and 68 nm for the 1201 bp one (data not shown). 
These high values might be due to the particle-substrate 
interaction, or particle-polymer or monomer-monomer 
excluded volume interactions which may swell the DNA. 


B. Refined extraction of the DNA persistence 
length using simulations of the HT-TPM setup 




FIG. 3: Simulated i?|| vs. Lp for DNA of (a) 2060 bp 
{N = 172) and (b) 1201 bp {N = 100). Triangles correspond 
to a free polymer-particle complex without excluded volume, 
full (resp. open) circles to a grafted polymer in the TPM 
geometry with a hard core particle and with (resp. without) 
intra-chain excluded volume interactions (see text). The solid 
line is the discrete WLC formula, Eq. (]^, the dashed one is 
the continuous WLC one, Eq. The dotted lines are fits by 
polynomials: (a) R|| = —10.7 x 10“®Lp-|-2.71 Lp-I-162.1 [nm], 
(b) Aji = -9.0 X lO-^L^ -h 1.955Lp -b 137.35 [nm]. 


To check these effects and the approximations used in 
Eqs. ( |3|4| ), we performed numerical simulations. DNA- 
particle conformations are computed numerically by ex¬ 
act sampling [23]. The labeled DNA polymer is gen¬ 
erated as a random walk of N steps, corresponding to 
the links of length 2a, with a bending energy by step, 
Llbend = —Kb COS 0, where Kb is the bending modulus and 
6 is the angle between successive steps. The number of 
steps N was chosen such that a = 6 bp, correspond¬ 
ing to the DNA diameter. The starting point is on the 
substrate and at each step self-intersecting trajectories 
(respectively trajectories intersecting the substrate) are 
eliminated to take into account intra-chain excluded vol¬ 
ume interactions (resp. repulsive interaction with the 
substrate). Hence the polymer is modeled by a chain of 
beads of (excluded) volume v = 47ra^/3 ~ 36 nm^, taken 
constant as a function of the ionic strength (see Fig. [^. 
The salt effects are therefore supposed to be completely 


taken into account in the bending modulus, Kb- The last 
step of length a + Rp has a uniformly random orienta¬ 
tion. The persistence length value is related to Kb by 
Lp = 2aPKb, where (3 = {kBT)~^ is the inverse of the 
thermal energy m- The two-dimensional projection of 
the particle position, R| |, was measured and the ampli¬ 
tude of motion, defined as i?|| = y'^(R^), is averaged 

over several millions independent trajectories. Since at 
this level of coarse-graining, electrostatic interactions are 
not included, we varied Kb by hand such that Lp spans 
the range 36 to 70 nm, in order to reflect the stiffening 
due to the decrease of the ionic strength in the solution. 

Figure shows the simulated i?|| as a function of Lp 
(full circles). We observe an increase of i?|| from 240 
to 300 nm for the 2060 bp long DNA, and from 195 to 
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230 nm for the 1201 bp one. Both ranges of Rn val¬ 
ues contain the corresponding experimental observations, 
which thus indicates that we explored the good range of 
persistence length values. 

For the purpose of comparison are also plotted the 
end-to-end distance of a free polymer-particle complex 
without any excluded volume interactions and without 
wall (triangles), and the polymer-particle plus substrate 
without excluded volume (open circles). The solid curves 
correspond to the discrete WLC result for the end-to-end 
distance, using Eq. (|^ (to include the particle contribu¬ 
tion) and without any excluded volume and wall ^5] : 


(Rf|) = ^[a2iVfF^(u(«6))+i?2] 


where 


Wn{x) 


1 -I- X 2x 1 — 

1 — X N {1 — xY 


(5) 

( 6 ) 


and v{kii) = coth(Kf,) — 1/Kb- Clearly the solid curve per¬ 
fectly matches the simulation results, as expected. The 
dashed one corresponds to the continuous WLC (no ex¬ 
cluded volume and no wall), Eq. Q, which gives slightly 
larger end-to-end distances. 

One observes that the presence of the particle which 
interacts both with the substrate and the chain induces a 
non-constant shift to higher values of i?|| (from triangles 
to full circles). Since the intra-chain excluded volume 
swells the polymer by less than 2 nm, especially for small 
values of Lp (more flexible chains), the main difference 
comes from the substrate-particle interactions. This is 
the reason why the extraction of Lp using Eqs. ffl or 
equivalently Eq. ([^ overestimates Lp of about 20 nm. 

To obtain precise values of Lp from the experiments, 
we thus fitted the simulation data R\\{Lp) by a quadratic 
polynomial law (see Fig.[^, which in turn allows us to ac¬ 
curately determine the experimental Lp from the experi¬ 
mental values of i?||. The persistence length is then plot¬ 
ted as a function of I~^ in Fig.j^for the two DNA lengths 
and the two types of counter-ions (figures are given in 
Tables 1 and 2 of the Supporting Information). Other 
available data found in the literature are also shown. 

With Na+ counter-ions, Lp values are in the same 
range both for L = 2060 bp and 1201 bp, which tends to 
confirm that the persistence length extracted with this 
procedure is almost independent of the DNA length, as 
expected. It increases monotonically from roughly 35 nm 
for high ionic strength (/ ~ 3 mol/L) to 54 nm for low 
one (/ = 10 mmol/L), which corresponds to an increase 
of more than 50%. Near physiological salt conditions, 
around 150 mmol/L, we find Lp ~ 43 nm. Moreover 
using this plot representation, the data show a clear con¬ 
cave shape. 

With Mg^+ counter-ions, Lp is greatly reduced which 
is a signature of the role of the ion valency z. This has 
already been observed in previous experiments [HI?]. The 
Lp values, between 35 and 50 nm, are slightly different 
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FIG. 4: DNA persistence length, Lp, vs. the inverse of the 
ionic strength, I~^, extracted from the HT-TPM data (full 
circles for L — 2060 bp and open ones for L = 1201 bp), and 
with available data (Refs, are given in the legend), (a) Case of 
buffers with monovalent Na^ salt counter-ions, and (b) with 
added divalent Mg^"*" counter-ions. 


for the two DNAs, Lp being larger by almost 5 nm at low 
/ for the longest DNA. Moreover, we observe an abrupt 
decrease of Lp between the case of no divalent ions and 
the previous I~^ value, corresponding to the addition of 
0.15 mmol/L of Mg^+. At higher /, the increase of Lp is 
almost linear in I~^. 

In the following we shall try to fit the so-obtained DNA 
persistence lengths using the available theories found in 
the literature. Before this, we compare our experimental 
values to the ones found in the literature. 


IV. COMPARISON WITH OTHER 
EXPERIMENTAL DATA 

A. Influence of Na^ monovalent ions 


We start the comparison with other experimental stud¬ 
ies by focusing on the measurements of Lp for values of / 
close to the physiological salt conditions.The behaviour 
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of Lp under the action of ionic strength is then discussed. 

1. Persistence length close to the physiological salt 
concentration 

An inventory of the values of Lp measured near 
the physiological salt conditions, i.e. for / G 
[100, 200] mmol/L with only monovalent Na'*' counter¬ 
ions, is presented in Table The mean value of these 
Lp values found in Refs. 0 [THU HZ] is 48 ± 6 nm 
which is in good agreement with our interpolated value 
Lp(150 mmol/L) = 43 ± 3 nm. In fact, the measured Lp 
varies widely with the experimental techniques in use, 
the studied DNA which lengths varies from 6646 bp to 
50 kbp, as well as with the theoretical and analytical tools 
used to extract Lp. In addition, the variability of these Lp 
values might come from the difficulty to perfectly control 
the presence of divalent ions such as Mg^+, the presence 
of which can have a dramatic effect even at low concen¬ 
trations (on the order of mmol/L), as mentioned before. 
Finally, it can be noted that the commonly accepted 
value of Lp = 50 nm for a “random” DNA sequence 
at physiological salt conditions, i.e. / ~ 150 mmol/L 
{I~^ ~ 6.7 L/mol) with only monovalent Na'*" counter¬ 
ions [26] slightly exceeds the experimentally derived ones. 

2. Variation of Lp on the whole I range 

In order to compare the global behavior of Lp{I) that 
we measured with the previously published results, we 
superimposed all the results in Fig. where Lp values 
are plotted as a function of I~^. As mentioned in the 
Introduction, Savelyev [3| has recently reviewed all the 
available experimental data and showed that they could 
be divided into two groups based on the difference in 
Lp behaviours observed at high ionic strength, 0.11 < 
/ < 3 mol/L (see Fig. [^. We will keep this division to 
compare our results to those obtained by the first group 
of experimental data [DEHU that shows a slow decrease 
of Lp with increasing /, then to those obtained by the 
second one laiHHIH that found a significant one. 

In the first set of experiments, the authors observed a 
saturation of Lp at high I. More precisely for I exceeding 
150 mmol/L, Lp decreased by only 10% only. Moreover, 
at low /, these publications show an increase in Lp of 
about 10% similarly as we do observe on our HT-TPM 
measurements. It can be noted that Harrington [Tj found 
a larger increase of Lp in the low / range than all the 
other results. 

This first group of experimental works gathers re¬ 
sults obtained by FB of Harrington [Tj , MB of Maret et 
al. [Bj [BS] , FOT of Baumann et al. [7] , and LD of Rizzo 
et al. [5]. The FB, MB or LD methods may be prone to 
perturbative Joule heating and bulk electrophoresis ef¬ 
fect; to minimize them, the ionic strength was therefore 
kept low. As FOT method induces a perturbation of the 


I 

[mmol/L] 

Lp 

[nm] 

Ldna 

[bp] 

Experimental 

method 

Ref. 

223.0 

66 ±3 

6646 

FB (25°C) 

m 

210.8 

46.5 

6646 

DLS (20°C) 

[SI 

210.8 

54.7 

6646 

DLS (20°C) 

m 

210.8 

40.6 

6646 

DLS (20°C) 

m 

204.1 

48 ±2 

39936 

DLS (20°C) 

m 

201.0 

40.6 ± 0.4 

39936 

FB (25°C) 

CTO] 

201.0 

46.8 ±0.4 

39936 

LD 

0 

186.2 

56 ±3 

48502 

FOT (25°C) 

m 

165.1 

40.5 ±0.4 

2060 

HT-TPM (25° C) 

this work 

165.1 

39.5 ±0.5 

1201 

HT-TPM (25° C) 

this work 

154.0 

50 ±5 

434, 587 

TEB (20°C) 

El 

123.0 

74 ±3 

6646 

FB (25°C) 

m 

110.1 

47.8 ±0.7 

2060 

HT-TPM (25° C) 

this work 

110.1 

46.2 ±0.8 

1201 

HT-TPM (25° C) 

this work 

103.1 

44.6 

6646 

DLS (20°C) 

m 

103.1 

44.6 

6646 

DLS (20°C) 

m 

103.1 

53 ±2 

39936 

DLS (20°C) 

m 

102.4 

55 ±2 

39936 

DLS (20°C) 

m 

101.1 

43 ± 1 

43 - 179 

TED (20°C) 

m 

101.0 

47.8 ± 0.4 

39936 

LD 

0 

101.0 

41.8 ±0.6 

39936 

FB (25°C) 

m 

93.4 

43 ± 1 

48502 

FOT (25°C) 

m 


TABLE I: Summary of Lp measured in the close-to- 
physiological salt conditions (with Na"*") found in the lit¬ 
erature. In order to compare the whole set of published 
data [BUTHSIIIZ], we rigorously computed the ionic strength 
used in these data by directly taking experimental values 
when available, or the values deduced by interpolation other¬ 
wise. 


sample structure in the high force regime, we only con¬ 
sider the Lp values obtained in the low force regime using 
the inextensible WLC model as a comparison to our HT- 
TPM results, where no force is applied. It is important 
to note that these studies were performed on DNA ex¬ 
ceeding 40 kbp long, and that only scarce measurements 
were performed at high ionic strength, in opposition to 
the second group of data. 

In this regime of high ionic strength, the second group 
did not show any plateau but rather a significant de¬ 
crease of Lp of about 25 to 30%. This observation is 
in good agreement with our measured decrease of 25% 
in the same I range. In the low ionic strength range, 
0.01 < I < 0.1 L/mol, this second group measures a but 
regular decrease of Lp of about 15% when / increases 
which is larger than the one we measured, about 10%. 
FB data m can be classified in this set of data, given 
the 25% increase of Lp in the high I range. However, 
variations in the low ionic strength range seems to be 
modest in comparison to the other publications of this 
set but closer to ours. 

This second group of experimental data is essentially 
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based on the DLS technique performed by Sobel et al. [T^ 
and Kam et al. (1982) [9]. Manning [8] and Post m 
proposed corrections to the extracted Lp values from the 
original data of Borochov et al. [3] . In these DLS exper¬ 
iments, the DNA hydrodynamic radius is deduced from 
the diffusion coefficient measurement. To infer Lp, the 
usual Gaussian polymer model was used. The number 
of Kuhn segments being large for the 6646 bp (resp. ~ 
40 kbp) DNA under study, A^ ~ 22 (resp. N ~ 1333), the 
swelling of the chain was induced as a result of excluded 
volume. Therefore, estimating precisely the excluded vol¬ 
ume is essential. It is nevertheless a challenging task. For 
instance. Manning’s [5] and Post’s m correction led to 
Lp values differing by 4 nm at / = 8 mmol/L, to ^ 1 nm 
at / = 1 mol/L. 

Note that we compare studies performed on long DNA, 
from 6646 bp to 40 kbp, which are in the flexible regime 
and thus much more sensitive to excluded volume effects 
than our HT-TPM experiments made on DNA of lengths 
L — 2060 or 1201 bp, which are in the semi-flexible 
regime. 

In this quantitative comparison, we only considered ex¬ 
periments perfomed on the same range of ionic strength 
induced by Na"*" ions. As a result. Refs. [H I33H53] that 
studied the effect of very low Na“''-induced /, as well as 
Refs. [SH [3S] that focused on the influence of K+ or Li+ 
ions, or Ref. [36] that studied the combined effect of mul¬ 
tivalent and monovalent ions, were discarded from our 
comparison. Note that single molecule study in Ref. m 
was not considered either due to the too few values of / 
explored. 


B. Influence of divalent ions 

The decrease of Lp, that is measured when I varies 
from 10 to 310 mmol/L, is about 35% whatever the ion 
valency and the DNA length (Fig. |^. Yet the shape of 
this decrease as a function of I is however completely 
different depending on the valency of the ions used. Di¬ 
valent ions appears to induce an almost linear decrease 
of Lp as a function of when monovalent ions cause a 
decrease with a concave shape. In addition, the absolute 
value for Lp is smaller by about 5 nm with divalent ions 
than with monovalent ions at low I. This observed trend 
is in agreement with the first quantitative observations 
at very low I by Hagerman |2], Elias and Eden [28], and 
then by others [lUj. In addition, Hagerman [3] observed 
on 434 and 587 bp DNA by TEB, the same abrupt fall of 
Lp at very low I with Mg^+ ions followed by a slower de¬ 
crease, as shown in Fig. [^ The magnitude of this initial 
decrease was however larger by around 30% in [5], 60% 
in [7] whereas [5] found a decrease similar to ours. 

Some other studies explored the influence of Mg^+ on 
DNA flexibility [27l jST] [33l [35| |39l [40| and showed a rise 
in DNA flexibility with the addition of Mg^'*' counter¬ 
ions in solution. Nevertheless, these studies only probed 
a few values of I, and some KCl was added to the solution 


preventing any quantitative comparison. 

Dietrich et al. [38| also used the TPM technique to 
monitor the effect of Mg^+ on a DNA fragment of 4882 bp 
and observed a large decrease of i?cxp|| in presence of di¬ 
valent ions. Yet Lp values extracted from these experi¬ 
mental data, much smaller than any other published ex¬ 
perimental data, cannot be quantitatively compared to 
our results due to several errors in the extraction pro¬ 
cedure. The persistence length was extracted assuming 
that the particle excursion was related to L using a simple 
Hooke law in the Gaussian regime, (R-^|) = ‘^LpL, thus 
seemingly forgetting the factor 3/2 due to dimensionality 
and not subtracting the particle radius [see Eq. ([^], but 
also ignoring excluded volume effects. Moreover, they 
presented a very large error bar of 10 nm, due to the 
small number of trajectories, ranging between 6 and 27. 

In the next section, we compare our experimental val¬ 
ues of the DNA persistence length to the various theories 
developed in the literature. 


V. COMPARISON BETWEEN EXPERIMENTS 
AND EXISTING THEORIES 

A. Odijk-Skolnick-Fixman approach at high ionic 
strength 


Several models have been proposed in the literature 
to explain the variation of the persistence length Lp of 
polyelectrolytes with the ionic strength I. When elec¬ 
trostatic interactions between mobile ions and the poly¬ 
electrolyte are taken into account at the Debye-Hiickel 
level (mean-field level and approximation of small val¬ 
ues of the electrostatic potential), it has been shown by 
Odjik [13] and Skolnick and Fixman [14] (OSF) that, us¬ 
ing a perturbative approach around an inhnitely stiff rod, 
the persistence length has two contributions 




(■B 

4A^k2 


(7) 


where is the bare persistence length (in the limit 
kA —> oo), and the second term is an electrostatic contri¬ 
bution to Lp, where A is the distance between elementary 
charges along the chain and 

K = (87r£B/)^/^ (8) 


is the Debye-Hiickel screening parameter. The Bjerrum 
length, is = e"^/{ATreksT), is equal to 0.715 nm in wa¬ 
ter at room temperature which yields n = 3.29-\/7 nm“^ 
(where I is in mol/L). Eq. Q is theoretically valid for 
polymer conformations close to the rod-like one, i.e. for 
£bL^ > A2. 

For a dsDNA, with two phosphate anions per bp, we 
have A = 0.17 nm, and « 50 nm so that the validity 
of Eq. Q is well verified. As a function of the ionic 
strength I, Eq. 0 can be rewritten as 


LPSF ^ 


0.559 


(9) 


I 




where / is in mol/L. Note that the numerical value of 
the coefficient depends on the model as shown by Fix- 
man m- In any case, Fixman affirms that the ex¬ 
ponent in I~^ is robust at large I (large k). Clearly 
the data shown in Fig. where the experimental Lp 
is plotted vs show a concave shape for high ionic 
strength, and therefore are not well fitted by the linear 
law, Eq. In Fig. is shown a linear fit (dashed 

lines), -I- C/I, for 0 < < 10 L/mol, which 

yields C = 1.21 nm-mol/L and = 35.5 nm for 
L = 2060 bp, and C = 1.52 nm-mol/L and = 31.7 nm 
for L — 1201 bp, which are almost 3 times larger than 
the one predicted by OSF. Since we do not have enough 
data at very high I with added Mg^+, we did not try to 
fit Eq. (|^ for the Mg^+ case. 

The data obtained by the first group II1I5H7] at high I 
were qualitatively in agreement with the OSE theory, in 
the sense that the variations of Lp are small at high I. 
However, no fit using Eq. 0 was done in these works. 


B. Manning charge renormalization at low ionic 
strength 



FIG. 5: Linear fits of the persistence length Lp measured 
by HT-TPM (already shown in Fig. vs. the inverse of 
the ionic strength . Dashed lines correspond to Eq. Q 
for large 7“^, 7//° -|- C/7, with 2 fitting parameters (see text 
for values), and solid lines to Eq. (Ill for low 7“^, with L® 
as fitting parameter. Better fits are obtained for Mg^^ by 
leaving B free (dotted lines, Lp = 38.5 nm, B = 0.121 for 
L = 2060 bp, and Lp = 36.5 nm, 73 = 0.1 for L = 1201 bp). 


The OSF theory fails to explain the variation of 
Lp{I~^) for the whole range of 7 studied. Indeed it does 
not reproduce the concave shape observed in Fig. A 
plausible explanation is that the Debye-Hiickel approx¬ 
imation of small values of the electrostatic potential is 
not valid for double-stranded (ds)DNA, which is highly 
charged. Indeed, Manning has shown in the 60’s that, 
in the limit of low 7, some counter-ions somehow “con¬ 
dense” close to the DNA due the large net charge of 
DNA [15]. This effect is a non-linear effect associated 
with the mean-field Poisson-Boltzmann equation close to 
charged cylinders. This phenomenon, known as the Man¬ 
ning condensation, tends to reduce the effective charge of 
the DNA. Manning proposed that, in the low salt limit, 
the effective linear charge density is a/A where a depends 
on the parameter 


u = 


zIb 


( 10 ) 


where z is the counterion valency: If u < 1, a = 1 and if 
M > 1 then a = l/u. For dsDNA, the Manning parameter 
is It = 4.21 for Na+ counter-ions (z = 1) and u = 8.41 for 
Mg^+ ones (z = 2). Thus for monovalent counter-ions, 
a = 0.24 and the effective charge is decreased by a factor 
of roughly 75%. For a mixture of counter-ions with differ¬ 
ent valencies, which is the case in our experiments with 
Mg^+ ions, the effective charge is a = A/{zIb) where 
z is the largest valency, i.e. only counter-ions with the 
largest valency (here divalent ones) condense along the 
DNA and no monovalent ions are condensed (unless di¬ 
valent counter-ions are depleted, which is not the case in 
most of the experiments and especially in ours) |43j . 

The OSF equation, Eq. ([^, has thus been modified 


according to 


rOSFM _ rO I 1 

P P 4z27bk2 


( 11 ) 


which changes the slope 73 of the linear relationship 
Lp = Lp + B/I from 0.559 [Eq. (j^j to 0.033 nm-mol/L, 
and should be valid at low ionic strength, i.e. for large 
This is the reason why the constant Lp is a pri¬ 
ori different from the one of Eq. (j^, L^. Note that 
adopting the Manning condensation, valid at low 7, to 
the OSF calculation of the persistence length, valid at 
high I is somewhat not consistent. Moreover Eq. (11) is 
not a proper asymptotic expansion in the limit kA —>■ 0 
since it leads to a diverging Lp. Of course, this limit can¬ 
not be reached in practice, since DNA counter-ions and 
ions resulting from water dissociation are always present 
even when no salt is added, which ensures that 1^0 
(even if it can be very small). Theoretically, the salt- 
free case corresponds to k~^ 3^ L, i.e. to ionic strengths 
7 ^ 10“® mol/L, which is far from being the case in our 
experiments. 

We fitted Eq. ( |Il| ) on our experimental Lp values at low 
salt for the Na+ case keeping as a free parameter. The 
results, shown in Fig.[^(solid lines), are quite satisfactory 
with Lp = 51.3 nm for L = 2060 bp and 48.8 nm for 
L = 1201 bp (letting 73 free leads to a slightly higher 
value of 73 = 0.038 nm-mol/L). However Eq. 0 does 
not fit the Lp values for the Mg2+ case, the slope being 
larger (we found 73 ~ 0.1 nm-mol/L) whereas Eq. ( [II| ) 
predicts a slope divided by z^ = 4, 73 = 0.008 nm-mol/L- 

The persistence lengths measured in [Millie] at low 
7 with Na'*' ions were fitted by Lp = L^ -|- 73/7. In [7] 
73 = 0.033 nm-mol/L was fixed to the Manning value. 
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and Lp was found to be 45 to 50 nm, which is com¬ 
parable to our results. Maret et al. fixed Lp = 
50 nm and obtained reasonable fits with 0.024 < B < 
0.041 nm-mol/L, which suggests a large error bar on 
the experimental values. Rizzo et al. also fitted 
their data for 3 < / < 1000 mmol/L. They obtained 
=46±1 nm in agreement with our value (see Fig. [^, 
and B = 0.043 nm-mol/L. Wenner et al. [35] used mea¬ 
sured Lp of ds DNA for various Na"*" concentrations by 
fitting force-extension curves at low forces. They ob¬ 
tained Lp = 46 nm and B = 0.037 nm-mol/L. All these 
values of B are in agreement with ours and the one pre¬ 
dicted by Eq. (Ill, 0.033 nm-mol/L. 

Tomic et al. m did dielectric spectroscopy experi¬ 
ments on semi-dilute DNA solutions with NaCl to inves¬ 
tigate the high frequency and low frequency relaxation 
modes vs. added salt concentration strength. In the 
high added salt limit (and relatively low DNA concen¬ 
tration) the length scale of the low frequency relaxation 
mode, Llf; can be interpreted as the DNA persistence 
length Lp. Their results are in qualitative agreement 
with the OSF-Manning theory, but with a coefficient 
B — 0.08 nm-mol/L larger than the Manning value, and 
smaller that the OSF one. Note that this discrepancy 
can be due to the fact that, in these experiments the 
total ionic strength is different from the added salt con¬ 
centration. Using the same experimental method, they 
investigated the effect of Mg^+ in solution [35] . The Llf 
was about 1.5 times shorter in Mg-DNA solution than in 
Na-DNA solution, also suggesting an increased screening 
with Mg^+. The behavior of Llf was again explained by 
the OSF-Manning theory but with a different value of 
the effective linear density. 


Later, Manning proposed to modify Eq. 0 by multi¬ 
plying the salt-dependent persistence length by a factor 
{2u—l)/u ~ 1.76 for 0 = 1 and 1.88 for z = 2 which gives 
a worse result for Na’*' and is not sufficient for Mg^+ in 
our case [49]. This correction is therefore not suitable. 


C. Mean-field non-linear corrections at 
intermediate ionic strengths 


In any case, Eq. (11) does not explain the concave 
shape shown in Fig. [7 at intermediate ionic strengths. To 
do so, one must have a renormalization factor of the DNA 
charge, a that depends on the ionic strength I. Such ap¬ 
proach has been developed by Netz and Orland m , in 
which the Poisson-Boltzmann equation is variationally 
approximated by a Debye-Hiickel one with a as a vari¬ 
ational parameter which renormalizes the electrostatic 
potential at the DNA surface. 

Here, we do the calculations by assuming that the 
DNA is not penetrable to ions, contrary to Ref. [50] . 
The dimensionless electrostatic Debye-Hiickel potential, 
(j) = zeP'tjj, for an electrolyte (valency z) of Debye-Hiickel 
constant k around a cylinder of radius R and surface 



FIG. 6: Renormalization charge parameter a vs. dimen¬ 
sionless screening parameter kR, solution of Eq. (151 for 
u = 4.11 (monovalent counter-ions in red) and u = 8.23 (di¬ 
valent counter-ions in blue). The black circles correspond to 
the experimental values studied in Section III. 


charge density a = (2ttAR) ^ empty of ions is: 


(/(r) = 


2u Ko{Kr) 
kR Ki{kR) 


for r > R 


2u Ko{kR) 

(hir) = -;-- lor r < R 

kR Ki{nR) 


( 12 ) 

(13) 


where Kq and Ki are the modified Bessel function of 
order 0 and 1 and u = zisjA is the Manning parameter. 

Following Netz and Orland [50], the full non-linear 
Poisson-Boltzmann equation is solved variationally by 
assuming that the solution is of the Debye-Hiickel form, 
a(j){r), where (/(r) is given in Eq. (12), and a, the fraction 
of the “free” counter-ions, is the variational parameter, 
and is solution of (see the Supporting Information): 


Airzisi^ — a) J drp{r)cj){r) 


J (irH(r)(/(r) {sinh[a(/(r)] — a^(r)} (14) 

where p(r) = a5{r — R) is the charge distribution. Note 
that the ionic exclusion factor D(r) was incorrectly put 
just in front of the sinh in [^. Using Eq. (121, Eq. (14) 
simplifies to 


2u(l - a)Ko{KR) = 

^ kR Ki_{KR)_ 


/ xKq{x) < sinh 
JkR I 


2ua Kq{x) 
kR Ki{kR) 


(15) 

dx 


The solution a{KR), where i? ~ 1 nm is the DNA radius, 
is plotted in Fig. for 2 = 1 (red) and z = 2 (blue). 
The renormalization factor a is a monotonous increasing 
function of kR with a^nR —>■ oo) —>■ 1 and a{KR —>■ 0) —>■ 
1/u. Hence it induces a concave shape to Lp defined as 


r MF _ r 0 




[a{KR)Y 


(16) 
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FIG._^ Same as Fig. where das hed lines are fits using 
Eq. (161 and solid lines using Eq. (19l. Parameter values 
are given in the text. 


One observes in Fig. [^that the variations of a{KR) are 
greater for z = 2 than for z = 1. Moreover for kR <C 1, a 
for z = 2 is smaller than for z = 1. These two features are 
in qualitative agreement with what is observed in Fig. 

To do a quantitative comparison of Eq. (161 with the 
experimental data in the whole range of ionic strengths, 
we used a polynomial interpolation fonction to fit Lp for 
monovalent ions (z = 1), and a power law for divalent 
ions (z = 2) ainR) ~ 0.423 {kR)^'^^‘^, shown in Fig. 
The fitting parameters are Lp and a prefactor in front 
of the second term of the rhs. of Eq. ([T^ (expected to 
be close to 1). Fits are shown in Fig. m as dashed lines. 
Clearly this approach leads to a slightly concave curve 
for Lp{I~^) for the two types of counter-ions. However, 
whereas fits of Mg^+ data (in blue) are reasonably good 
(LO = 35 .7 nm, prefactor equal to 1.8 for L = 2060 bp, 
and Lp = 35.2 nm, prefactor of 1.5 for L = 1201 bp), 
the fits of the Na+ data are not good for I > 0.1 mol/L 
(/“^ < 10 L/mol), the concavity being not enough pro¬ 
nounced {Lp = 44.2 nm, prefactor of 0.8 for L = 2060 bp, 
and Lp = 45.9 nm, prefactor of 0.9 for L = 1201 bp). 
Moreover, here again the non electrostatic contribution 
to the persistence length, Lp varies from the fit of Na"*" 
data to the fit of the Mg^+ ones. 


D. Beyond mean-field: ion-ion correlations and 
thermal fiuctuations 

Other approaches, that go beyond the mean-field ap¬ 
proximation by taking into account ion-ion correlations 
or/and thermal fluctuations, have been proposed [5TH55] . 

Nguyen et al. |53j considered ion-ion correlation in the 
strong coupling regime defined as F = 

1. The theoretical limit F —>■ oo corresponds to the freez¬ 
ing of the strongly coupled counter-ions into a Wigner 
crystal close to the DNA molecule. For F ^ 1, they 


found the following correction to the persistence length 
in the limit of zero ionic strength: 

d2 

Lp = LlL- —T(-0.83r -k 0.33r^/^ -k 0.87) (17) 

^ P 2zA 


which is independent of the ionic strength I. Therefore 
they propose that Lp is constant at vanishing I. For DNA 
at room temperature, one has F = 1.2 z^/^, and applying 
Eq. ( [TtI yields a correction of -1-0.65 nm for z = 1 (even 
if F ~ 1), and —2.20 nm for z = 2. Eq. ( [T7| ) therefore 
qualitatively explains the observed abrupt decrease of Lp 
of about 5 nm, when Mg^’*' ions are added at very low 
/ in the buffer: the monovalent Na’*' ions are replaced 
by strongly coupled divalent counter-ions, whose corre¬ 
lations decrease the global free energy of the condensed 
ions and therefore the DNA bending free energy. Hence 
Lp{I 0) decreases with z. Following this approach, it 
is thus consistent to choose two different asymptotic val¬ 
ues for Lp when / —>■ 0 for monovalent and divalent ions. 
This explains why these values where slightly different 
(of about 9 nm) in the preceding section. Note that this 
theory explains the constant shift at very low I but does 
not explain the change in the shape of Lp{I~^). 

Thermal fluctuations have been taken into account by 
Golestanian et al. [S2], who obtained a correction to the 
persistence length due to fluctuation-induced correlations 
between ions. Indeed, they correct the OSF-Manning 
formula Eq. (Ill at low kA, following 


L - L° \ 

P 4m2(kA)2 


f{KA,u) 


(18) 


where f{KA,u) = [1 — 2{u — 1) In(KA)] Ariel 
and Andelman m proposed a similar correction with 
/(kA, u) = u{2 — u) — {u— l)^/[Mln(tt;A)]. 

Both formulae do not apply to DNA for z = 1,2 
{u = 4.21 and 8.41). Golestanian’s formula yields a slope 
which is almost divided by 200 with a convex shape for 
< 10 L/mol, and Ariel’s one yields a decreasing func¬ 
tion of Lp{I~^) as increases for the whole I range. 
These two theories are therefore not consistent with the 
whole set of experimental data shown in the experimental 
results Section. 


E. Manning’s internal stretching force calculation 
for Ld 


The major issue in trying to fit the above theories for 
the whole range of ionic strengths is to find a fit which 
yields the concave shape observed for the experimental 
values. Eq. ([^ a priori valid for high I and Eq. (11) valid 
for low I cannot be reconciled since the constant value 
is clearly different in both cases, Lp > L^ and should 
therefore vary with I. 

Manning noted this discrepancy in 2006 [15] and pro¬ 
posed a new formula for the persistence length by tak¬ 
ing into account the internal electrostatic tension due to 
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the repulsion between charges along the polyelectrolyte. 
He adapted the calculation by Netz [54] for strongly 
stretched polyelectrolytes at the Debye-Hiickel level to 
the framework of the counterion condensation approach 
to obtain the persistence length of a polyelectrolyte as a 
function of kA and the persistence length of the so-called 
null isomer (the hypothetical structure of the polyelec¬ 
trolyte if the backbone charges are set to zero), L*: 


= Ca'-) 


2/3 ^ 4/3 


z‘^£b 

-l-ln(l-e-"^)] 


{2u - 1) 


kA e 
1 - 


(19) 


Eq. (191 fits very well our data for the Na’*' case on the 


whole ionic range with only one fitting parameter, L* 
(Fig. 0, The fitting values are L* = 6.0 nm (for L = 
2060 op) and 5.4 nm (for L = 1201 bp), close to the value 
of 7.4 nm fitted by Manning on Baumann’s experimental 
data [TK] . 

Savelyev |4] performed numerical simulations to inves¬ 
tigate the dependence of the persistence length of double- 
stranded DNA on solution with various ionic strength. 
A coarse-grained model of two-bead DNA chain with 
explicit mobile ions (Na"*" and Cl“ ions) [16] was de¬ 
signed to reproduce physical salt conditions from 10“^ 
to 0.1 mol/L (the water solvent is implicit). Their nu¬ 
merical results for Lp{I) are in semi-quantitative agree¬ 
ment with Eq. (19) for / > 0.1 mol/L. For lower / the 
agreement is better with the OSF theory, Eq. 0 (see 
Fig. 2 of Ref. [1], the fitting parameter values are not 
given). Moreover, Savelyev [4] compared previous exper¬ 
imental results to his simulations and found a qualitative 
agreement. 

Assuming that, according to the counter-ion conden¬ 
sation theory |4d] . all the condensed counterions are 
divalent, we use the same formula for the persistence 
lengths with Mg^+ counter-ions. It leads to poorer fits 
(Fig. [^ with very different values for L*: 14.1 nm (for 
L — OT60 bp) and 12.7 nm (for L = 1201 bp). The fact 
that L* varies, and increases, with the counter-ion va¬ 
lency is puzzling. Hence the dependence of Lp with z in 
Eq. (19) is not consistent with our experimental data. 


VI. INTERPOLATION FORMULA FOR THE 
WHOLE IONIC STRENGTH RANGE 


We propose the following interpolation formula to fit 
the four sets of data (2 DNA lengths, monovalent and 
divalent salts) on the whole / range: 


r — TOO 


r 0 r oo 

1 + {i/hY 


( 20 ) 


with 4 fitting parameters, Lp, Jq, and <5. The 

fits, shown in Fig. [^ are very good for the monova¬ 
lent Na"*" ion, with = 58.1 nm, = 33.8 nm, 
/q = 0.104 mol/L, and 5 = 0.931 for L = 2060 bp, and 



FIG. 8: DNA persistence length, Lp, vs. the inverse of 
the ionic strength, I~^, extracted from the HT-TPM data. 
The red (resp. blue) symbols correspond to buffers with 
sodium (resp. magnesium) counter-ions (top curves are for 
L = 2060 bp and bottom ones for L = 1201 bp). The solid 
lines are fits using Eq. (201 with parameter values are given 
in the text. 


Lp = 53.9 nm, L“ = 30.1 nm, Iq = 0.174 mol/L, and 
5 = 0.994 for L = 1201 bp. 

Several comments are in order. First, as expected, 
the asymptotic values, L^ and Lp, are, at about 2 nm, 
the same as the ones extracted from the linear OSF and 
OSF-Manning fits shown in Fig. Next, the values 
of the crossover ionic strength, Jg, is on the order of 
0.1 mol/L which corresponds to a Debye screening length 
Kq ^ ~ 1 nm, i.e. on the order of the DNA radius, R. 
It thus suggests that the concave shape, which is more 
pronounced for I ^ Iq comes from non-linear Poisson- 
Boltzmann effects as illustrated in Section C. Moreover, 
the effective power law for / ~ /g is found by doing a 
logarithmic expansion of Eq. (20) around /g: 


/ tO _ Too 

ln(Lp-L~)cslnMh^ 


-^n(/) (21) 


which yields an exponent —6/2 ~ 0.5 which is a good 
approximation for 0.05 < / < 0.5 mol/L (error less than 
1 nm). Finally, at low ionic strength, J <C /g, Eq. (20) 


varies slightly and the curve looks like a linear law as 
a function of with a small slope, as suggested by 
Eq. At large Ionic strength, J ;g> Jg, Eq. ( pO] ) yields 
Lp ~ L/° -I- (Lp — L^){Io/IY which is equivalent to 
Eq. (0 for i5 = 1 but with a slightly larger slope of 2 
to 4 nm-L/mol. 

For divalent Mg^+ ions, the fits are also good but with 
different values for 6 and /g as compared to the Na+ only 
case. The parameter values L/° are kept equal to the Na+ 
case, and the values of Lp are comparable: Lp = 59.2 nm, 
Ig = 0.017 mol/L, and <5 = 0.830 for L = 2060 bp, and 
Lp = 51.7 nm, /g = 0.045 mol/L, and 6 = 0.546 for L = 
1201 bp. As expected the concavity is less pronounced 
and shifted to lower values of I, close to /g. The formula. 
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Eq. (20), can be useful for experimentalists to interpolate 
values of Lp on the whole range of /. 


VII. CONCLUSION 


A. Summary 


Using the High-Throughput Tethered Particle Motion 
setup, we measured the impact of the ionic strength 
on DNA conformation for two DNA of lengths 2060 bp 
and 1201 bp. To this end, we investigated a large 
and homogeneously distributed range of ionic strengths, 
I G [0.01, 3] mol/L, by adding salt to the buffer with 
monovalent Na"*" or divalent Mg^’*' counterions. Exper¬ 
imental drift and biases due to the finite exposure time 
of detector were were corrected. To extract properly the 
DNA persistence length, Lp, from the HT-TPM ampli¬ 
tude of motion, iZoxp||) numerical exact sampling simu¬ 
lations (without explicit mobile ions and solvent) were 
performed. Both the DNA excluded volume and the par¬ 
ticle one’s were taken into account. These simulations 
allowed us to obtain the experimental Lp as a function 
of / with a good accuracy of about 4% (Fig. [^. When 
Lp is plotted as a function of the overall trend is a 
monotonous increasing function with, for the Na'*" case, 
a concave shape, and, for the Mg^+ one, an almost linear 
one (except at very low concentrations of Mg^+). 

Our results are compared to other results found in the 
literature. A quantitative comparison is difficult, since 
the Lp values fluctuate appreciably depending on the ex¬ 
perimental setup and the method of extraction of Lp. 
Hence, for instance at I ~ 150 mmol/L, Lp lies between 
40 and 74 nm (Table |T|. The available experimental val¬ 
ues can, as proposed by Savelyev et al. [1], be divided 
into two sets of data. The global behavior of our mea¬ 
sured Lp with I is not coherent neither with the first set 
nor with the second set of experiments. On one hand, 
our Lp values appear to be in agreement with the slow 
increase (of 10%) of Lp observed by the first group at low 
I. On the other hand, at high /, our Lp values show a 
significant 25% variation, in perfect agreement with the 
experiments of the second group. 

Our experimental Lp values follow a linear OSF law 
in only for a very small range of ionic strength at 
high I with a different prefactor than predicted by OSF 
in Eq. Q. 

For monovalent Na"*" counter-ions, Lp varies linearly 
with at very low I, according to the OSF equa¬ 
tion using the Manning counter-ion condensation theory, 
Eq. (11). The whole I range is furthermore well fitted 
by Eq. (19), which takes into account both the DNA in¬ 


ternal stretching due to phosphate ions of the backbone 
and the counter-ion condensation around the DNA. 

For divalent Mg^"*" counter-ions, however, neither 
Eq. (11) nor Eq. (19) (with the same fitting parame¬ 
ter as for Na'*') fit well the data. It suggests that these 
theories do not reproduce well the observed dependence 


with the valency z. Using a variational approach taking 
into account both non-linear Poisson-Boltzmann effects 
and screening by mobile ions, we propose a reasonable 
fit both for the Na"*" and the Mg^+ cases, but only for 
/ < 0.1 mol/L and with two different values Lp at van¬ 
ishing ionic strength. This marked decrease of Lp at very 
low I, when a very small amount of Mg^+ ions is added, is 
semi-quantitatively explained by a theory which consider 
ion-ion correlations, Eq. (0, for large z. 

In order to both interpolate Lp values in-between the 
ones effectively measured, and compare to future exper¬ 
imental data, we proposed an empirical formula which 
fits both the monovalent and divalent cases. 


B. Concluding remarks 

The large scattering of the available experimental Lp 
values observed in Figs. |^a) and (b) may be due to 
the different experimental setups but also to the vari¬ 
ous buffers used in these experiments. Indeed, we have 
shown that the presence of traces of divalent ions in the 
buffer can decrease substantially the Lp value at a given 
ionic strength. They are often present in buffers in order 
to maintain the fixed pH. 

An illustration of the extreme sensitivity of persistence 
length values to the experimental method and the model 
employed to extract the results was shown by Mielke et 
al. Hg. Brownian dynamics simulations on a double- 
stranded DNA in a bulk environment were performed 
at different salt concentrations. Two different strategies 
were employed to calculate Lp from the simulation re¬ 
sults. One used the expression Eq. ([^ of the WLC model, 
and the other an approximation proposed by Hager- 
man [2] of the rotational diffusion coefficients in order to 
directly connect to the experimental results of Ref. . 
At low concentrations, depending on the method used, 
Lp{I) splits in two distinct behaviours. Values from the 
rotational diffusion coefficients were more than 30 times 
larger than the WLC ones and Hagerman’s values. This 
result highlights the role of the chosen model to extract 
Lp, and the rough approximation used in earlier models 
to extract it from DLS measurements. 

Furthermore, it has been recently shown that HT-TPM 
can detect the effect of the DNA sequence, in particular 
the presence of A-tracts, on the DNA conformation and 
has been interpreted as a modification of the DNA spon¬ 
taneous curvature |19j . Preliminary results also show 
that, for a given DNA length but two different sequences, 
the persistence length varies. It is thus tempting to sug¬ 
gest that the bare, non-electrostatic, contribution to Lp 
can also be sequence-dependent and be another explana¬ 
tion for this data scattering. 

Finally, many experiments study either the influence 
of ions with higher valency, for instance trivalent ions 
such as spermidine which also has a strong effect on Lp 
at millimolar concentrations [3 135], and/or the role of 
multivalent ions on the DNA melting temperature [56j . 
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It would be interesting to pursue such a quantitative 
study of the DNA conformation for such trivalent ions, 
and thus to study the interplay between screening 
effects, condensation [57], and denaturation |58j . 

On the theoretical side, a complete theory which ex¬ 
plains the variations of Lp both as a function of / and the 
counter-ion valency z is still lacking. The approach by 
Manning m is appealing since it fits very well the exper¬ 
imental data for monovalent counter-ions with only one 
fitting parameter. However, the precise treatment of ion- 
ion correlations should be taken into account to extend 
such a theory to counter-ions with higher valencies. Note 
that these approaches do not consider dielectric exclusion 
close to a low dielectric molecule such as DNA [sniisni or 
van der Waals interactions which are also modified when 
the ionic strength is varied. 

In this paper we assumed that, according to the WLC 
model, the tangent-tangent correlation function is a sim¬ 
ple exponential, Eq. Q, which therefore leads to a single 
correlation length, Lp. However Barrat and Joanny show 
that, by taking into account the polymer bending fluc¬ 
tuations, the persistence length is scale-dependent m- 


A more appropriate choice would be a double exponen¬ 
tial where (t(s) • t(0)) ~ 1 — s/L'^ at small length scales 
and (t(s) • t(0)) = exp (—at large length scales, 
where the crossover depends on the ionic strength. It has 
been shown [62| that A®is given by Eq. Q and is 
the bare persistence length. Another approach proposed 
that the latter is also salt-dependent [B5] . Of course such 
a model is more difficult to apply to the experimental 
measure of i?DNA only, with a rather tricky extraction 
of two different correlation lengths, but it might provide 
a relevant framework to explain the overall observed be¬ 
haviours of Lp{T). 
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